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Abstract 

We derive a long distance effective action for space-time coordinates from a IIB matrix 
model. It provides us an effective tool to study the structures of space-time. We prove 
the finiteness of the theory for finite to all orders of the perturbation theory. Space- 
time is shown to be inseparable and its dimensionality is dynamically determined. The 
IIB matrix model contains a mechanism to ensure the vanishing cosmological constant 
which does not rely on the manifest supersymmetry. We discuss possible mechanisms 
to obtain realistic dimensionality and gauge groups from the IIB matrix model. 
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1 Introduction 



A large reduced model has been proposed as a nonperturbative formulation of type IIB 
superstring theory IQ 0. It is defined by the following action: 

^ = '^2Tr{\[A„ A,] [A^, AH + l^r>^[A„ ^]), (1.1) 

here ip is a ten dimensional Majorana-Weyl spinor field, and and are N x N Hermitian 
matrices. It is formulated in a manifestly covariant way which we believe is a definite ad- 
vantage over the light-cone formulation to study the nonperturbative issues of superstring 
theory. In fact we can in principle predict the dimensionality of space-time, the gauge group 
and the matter contents by solving this model. In this paper we report the results of our 
first such efforts toward that goal. 

This action can be related to the Green-Schwarz action of superstring by using the 
semiclassical correspondence in the large limit: 

Tr ^ N J d^a^g. (1.2) 

In fact eq.( p..l|) reduces to the Green-Schwarz action in the Schild gaugef^: 

^schiid = / rfV[^a(^{X^X^}2 - '-^T^{X^,ij})+pfg]. (1.3) 

Through this correspondence, the eigenvalues of matrices are identified with the space- 
time coordinates X^(cr). The Af = 2 supersymmetry manifests itself in S'schiid as 

(5(i)X^ = te\r^^, (1.4) 

and 

= e,, 

5(2)X^ = 0. (1.5) 
The M = 2 supersymmetry (|1.4|) and ( |1.5|) is directly translated into the symmetry of S as 

5(1) = leiT^i,, (1.6) 



and 

5^^^A^ = 0. (1.7) 

The bosonic part of the action vanishes for the commuting matrices = x^^*-' where i 
and j are color indices. These are the generic classical vacuum configurations. Namely the 
distributions of the eigenvalues determine the extent and the dimensionality of space-time. 
Hence the structure of space-time is dynamically determined by the theory. As we will show 
in this paper, space-time exits as a single bunch and no single eigenvalue can escape from 
the rest. We may consider n x n submatrices of our original N x N matrices where n is 
assumed to be still very large. More precisely, we consider a block diagonal background 
where the first blocks are n dimensional and the rest is a diagonal vacuum configuration. 
Such a background may represent a D-string (or D-objects) which occupies a certain region 
of space-time. 

Then the effective action which describes the n x n submatrices can be obtained by 
integrating the rest of the degrees of freedom. It has been shown that the term proportional 
to f3 in eq.( p..3|) arises if the vacuum energy of order exists in this model[|l|. In fact we will 
show in this paper that it is indeed the case. The effective action may include all possible 
terms which are consistent with the symmetry. The term proportional to (3 is the lowest 
dimensional term which is consistent with the symmetry and we expect that the higher 
dimensional terms are irrelevant for large n. In this way we can understand why this model 
resembles the effective action for D-objects 0. However we emphasize that we regard the 
action (|1 . 1| ) as fundamental and it should be distinguished from the effective action for 
D-instantons. 

It is straightforward to construct the configurations which represent arbitrary numbers of 
D-objects in space-time in an analogous way. Although they are constrained to be inside the 
extension of space-time as we see in this paper, otherwise they may be located anywhere. As 
is shown in the long distance interactions of D-objects are in accord with string theory. 
Thus it must be clear that the IIB matrix model is definitely not the first quantized theory 
of a D-string but the full second quantized theory. Arbitrary numbers of D-objects exist as 
local excitations of space-time in this formulation which is certainly not the case if we deal 
with the effective theory of a fixed number of D-objects. 

It has been also proposed that the Wilson loops are the creation and annihilation oper- 
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ators for strings. We consider the following regularized Wilson loops[0: 



w{C) 



Tr[v{C)], 



M 



v{C) 



n=l 



n 



(1.8) 



Here are momentum densities distributed along a loop C, and we have also introduced 
the fermionic sources A„ . e in the argument of the exponential is a cutoff factor. As it 
has been shown, e can be regarded as a lattice spacing of the worldsheet. We therefore call 
e — s> limit the continuum limit. This model has been further investigated through the loop 
equations which govern the Wilson loops. Such investigations have shown that the theory 
has no infrared divergences and the string perturbation theory can be reproduced in the 
double scaling limit. 

It is possible to see the infrared finiteness of the theory in a more direct way. We may 
expand the fields around the vacuum configurations = x^^*-'. While the off-diagonal 
elements of the matrices possess the nonvanishing propagators, the diagonal components 
do not possess the propagators and may be called the zeromodes. We denote the bosonic 
zeromodes as and the fermionic zeromodes as where i denotes the color indices. 

The perturbative expansion around this background at first sight appears to resemble 
that of the original large reduced models[0|. If so, the perturbative expansion around the 
classical vacuum is identical to the ten dimensional super Yang-Mills theory in the large 
limit. The ultraviolet divergences of ten dimensional super Yang-Mills theory appear as the 
infrared divergences in this model and hence the theory appears to be ill-defined. However 
the existence of the fermionic zeromodes turns out to modify such a picture completely as 
will be shown in this paper. 

The possibly dangerous configurations are such that the eigenvalues are widely sep- 
arated compared to the scale set by g"^. It is clear that the mass scale of the off-diagonal 
elements and that of the zeromodes are also widely separated in such a region. Therefore we 
can adopt a Born-Oppenheimer type approximation and we can integrate the off-diagonal 
elements first. Here the expansion parameter is g"^ over the fourth power of the average 
distance between the eigenvalues . Therefore the one loop approximation gives the accurate 
result for widely separated eigenvalues. In this way we can obtain the effective Lagrangian 
for the zeromodes. 
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We can subsequently perform the integration over the fermionic zeromodes. The resulting 
effective action for the bosonic zeromodes turns out to be completely infrared finite for finite 
N. Since the theory is manifestly ultraviolet finite, we have established the remarkable fact 
that this model is a finite theory which is free from both short distance and long distance 
divergences. 

The effective action provides us a useful tool to study the structures of space-time. In 
section three, we study it to determine the distributions of space-time coordinates. We 
believe that we have entered a new paradigm in string theory where we are beginning to 
determine the structure of space-time dynamically. The readily recognizable structure which 
emerges from our investigation is in fact a four dimensional object albeit it is a fractal. We 
find that space-time points form branched polymers in ten dimensions within the simplest 
approximation. Obviously, we need some mechanism to flatten it into four dimensions to 
describe our space-time. Although our investigation to find it in the IIB matrix model is 
still in progress, we explicitly propose possible mechanisms. 

We remark that related models to ours have been proposed and studied in |^~pO[. We 



also note a deep connection between our approach and the noncommutative geometry [21|. 
This paper consists of four sections. In the first introductory section, we have reviewed the 
conceptual framework of our approach. In section two, we derive the long distance effective 
action for the zeromodes (super coordinates). In section three, we study the effective action 
to understand the possible structure of space-time. Section four is devoted to the conclusions 
and the discussions. One of the crucial issues in our approach is how to take the large N 
limit. We address this question in the concluding section. 

2 Effective theory for diagonal elements 

In this section, we discuss an effective theory for the diagonal elements of the IIB matrix 
model. As we have explained in the introduction, the theory classically possesses the huge 
moduli and the diagonal elements of may assume any values. However this classical 
degeneracy is lifted quantum mechanically and there remains no moduli in the theory as we 
shall show in this section. Let us consider the expansion around the most generic classical 
moduli where the gauge group SU{N) is completely broken down to U{1)^^^. Then the 
diagonal elements of and ip appear as the zeromodes while the off-diagonal elements 
become massive. So we may integrate out the massive modes first and obtain the effective 
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action for the diagonal elements. 

We thus decompose into diagonal part and off-diagonal part A^. We also decom- 
pose ip into ^ and ip: 



A., 



X,, 



\ 



( e 



\ 



(2.1) 



where and satisfy the constraints Z)i=i 2;)^ = and '}2i=i Ca = O5 respectively, since the 
gauge group is SU{N). We then integrate out the off-diagonal parts A^ and and obtain 
the effective action for supercoordinates of space-time 5'eff[X, The effective action for 
space-time coordinates S'efrfX] can be obtained by further integrating out ^: 



/ 



dAdipe 



-S[A,i,] 



dXd^e^^^«^^'^^ 



(2.2) 



where dX and d^ stand for Hili^ dx^^ and Hili^ d^l^, respectively. 
2.1 Perturbative evaluation of Sqs[X,^] 

In this subsection, we integrate out the off-diagonal elements A^ and ip by perturbative 
expansion in g'^. As we will see below, this expansion is valid when all of the diagonal 
elements are widely separated from one another. Since X^ and ^ do not have propagators, 
we treat them as collective coordinates. 

The original action ( |1.1| ) can be expanded as follows. 



S 



So 



S2 + 'S'int 
1 



2g 



29' 



Jri-[X,,A,][X^,A^] + [X,,A,][X\A^] 

-i^T^[X„i,] - [lA,]T'^i,-i,T'^[A„^]), 
;Tr{-2[X^,A,][A^,A^] - 



-^r^[4,^]). 



(2.3) 



(2.4) 



(2.5) 
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We can observe that ^ has no propagator since the quadratic term of ^ is absent. As for the 
gauge fixing, we adopt the following covariant gauge: | 

^F.P. = -^Tr([X^,6][A^c]), (2.6) 

where b and c are the Faddeev-Popov ghost fields. 

In terms of the components, S2 + 5'g.f. can be written as 

Hc - e)rv'''A^'* + tp'^'^ic - e)^ii)- (2.7) 

The first and the second terms are the kinetic terms for A and respectively, while the last 
two terms are AipC, vertices. The basic building blocks of the Feynman rules are: 




<ri^''>^ =--— (2.9) 

■ (x*-x-')^r^ 



vertex = ^^^f^ ^ = r^(f - f ) (2-10) 

As A and ip are matrices, the propagators are denoted by double lines with indices i and j. 
The interaction terms in eq. ( |2.5|) are denoted by the following vertices: 



boson three-point vertex = ■ (2.11 




boson four-point vertex = / JT-n'C (2.12) 



k n j 

boson fermion vertex = (2.13) 
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Actually there is subtlety for the ghost zeromode in the covariant gauge. However the same one loop 
effective action can be obtained without such subtlety in the light cone gauge where A-^- is diagonalized. 
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One can see from the above Feynman rules that the expansion parameter after the ^ 
integration is g'^ over the fourth power of the average distances between the space-time 
coordinates. Therefore, the perturbative expansion is vahd when all of the diagonal elements 
are widely separated from one another. In particular, one-loop approximation is valid in these 
regions. 

A novel feature of the above Feynman rules compared with those of the gauge theory in 
the large N limit is the presence of ^ insertion vertices ( |2.10| ). If we were to set ^ = 0, the 
contributions from A and ip along from the ghosts cancel at the one loop level: 

dAd^dhdc e-(^2+s,.t.+Sp.p.) _ Y{{x' - a;^-)2(-io+2+8) = i. (2.14) 

Thus the diagonal elements would take any value and the classical moduli would remain at 
least perturbatively. The theory is identical to the large limit of the super Yang-Mills 
theory if we identify X with momenta^]. Therefore the major difference between the IIB 
matrix model and super Yang-Mills theory is the presence of dynamical zeromodes X and ^ 
in the former. Although the zeromodes are order A^ quantity, we cannot ignore them even 
in the large A^ limit as we shall show below. 

Eq. ( |2.7|) shows that the one loop effective action can be expressed by a super determi- 
nant. It turns out to be most useful to integrate out ip first. Such an integration gives rise 
to the following quadratic term of A^: 

i- _ ei)r'^-(e - e)|^f^ir^i^' (2-15) 

9 i<j [X -X-'} 

where T^"'^ = ^Fl^F^T'^l. Here the indices within the bracket [ ] are totally antisymmetrized. 
In this way we obtain the following one-loop effective action for the zeromodes: 

' dAd^dbdc e-(^2+5,.,.+5p.p.) _ -Q (iet^,{r]>"' + S^^.^)'^ 

i<j 

^ e-^cV-^[^'«], (2.16) 
where 

^cj) = - e^)F^"^(e - (2.17) 



The effective action can be expanded as 



q2 04 06 08 

-!:*.•(%+¥+%+%). (2.18) 



- 4 6 



Here the symbol tr in the lower case stands for the trace for Lorentz indices. Note that 
the expansion terminates with the eighth power of since ^ has 16 spinor components 
and S'f^^j^ is bilinear in ^. Note also that only even powered terms remain, since S^^j-^ is anti- 
symmetric under the interchange of fi and u. Furthermore the quadratic and sextic terms 
vanish identically as is shown in the appendix 0. Therefore we obtain 

c4 c8 

^eV°°'^[x,e] = -E^K^ + -f )• (2.19) 

i<j 



We remark that eq. ( |2.19| ) is gauge independent. It is because the longitudinal part of the 



bosonic propagator which is sensitive to the gauge parameter has no contribution since 

(< - ^i)s^^ = i^' - e^)r''-'^(e - = 0. (2.20) 



It is also possible to construct similar reduced models for super Yang Mills theory in 
D=3,4 and 6 by the large N reduction procedure. Then the calculation we have explained 
in this section is also applicable to these cases and we find 

r ^ for D=3 

Sls'°"''[X,^] = \-^i<jtr{^) ^ forD=4 . (2.21) 
i-Ei<,tK% + %) forD=6 

Before we proceed to perform C, integration and to obtain the effective action for X, we 
discuss Af = 2 supersymmetry of the effective action in the next subsection. 

2.2 A/" =2 SUSY 



In this subsection, we show that the effective action ( |2.19D has the following M = 2 super- 
symmetry: 

\ =0 

r 5(2)a;j^ = 

These symmetries are the remnants of the symmetry in the original theory: 



(2.22) 
(2.23) 



j 5(2)4 =0 

U(2V =62 ■ 



(2.24) 
(2.25) 
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First let us decompose the variables into the diagonal and off-diagonal parts in the original 
transformations ( p.24|) and (|2.25| ). We may consider and ip as the quantities of order 
h. Then one can expand these transformations in terms of h. However the action is invari- 



ant under these transformations at each order of h. Thus 5*2 
transformations which are linear in and ip, 



is invariant under the 



<5(i)iij- 





iix 










(2.26) 



After the integration over A and the remaining effective action ( |2.19| ) shall have the 
symmetries (O^) and (^1231). 



One can also show the invariance of (gl9D under the transformations (^^221) and (^J2^ ) 
through explicit calculations. Since S'^^--^ contains only through the combination as ^* — 
one can see that the invariance under ( p.23| ) is satisfied rather trivially. As for (|2.22|) , it 
introduces additional ^\ tr{S(ij)Y is invariant under ( p. 22 ), since it already contains 16 ^'s. 
Some calculations are required to exhibit the invariance of the term tr[S(ij)Y: 



1 1 



(2.27) 



where JJ^^'^p'' 



(eir/^Oler^'^^Oler/^Oler/Oler/^O is a tensor which is totally sym- 
metric and traceless with respect to the indices vXpa as is explained in the appendix A. Here 



and ^ denote (x* 



and ^' — respectively. We may further decompose f/^^^Vo" jj^to 



the irreducible components as 



sym. traceless 



symmetrization of U''"^'"' - — (^^"f/^"^'"' + 9 terms). (2.28) 



And we can show that t/sym.traccless 



by decomposing SO(IO) D S0(2) xS0(8) where 

s% = 0. 



{ha,ta). It is then sufficient to demonstrate that f^syi^it^!eiess 



Therefore the following cancellation implies the invariance of tr{S{,ij^Y 



[X 



(2.29) 



This completes the proof of the invariance of ( |2.19| ) under ( |2.22| ) and ( ^^.231 ). 
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2.3 Dynamics of X at long distances 

In this subsection, we study the integration procedure of ^ to obtain the effective theory for 
the space-time coordinates X: 

J (iXd^e-^^«[^'^l = J c/Xe-^=ff[^l (2.30) 

We perform ^ integrations exphcitly at the one loop level in this subsection and explain that 
it can be represented graphically in terms of trees which connects the space-time points. Al- 
though the vacuum energy of order A^^ has been canceled between the bosonic and fermionic 
degrees of freedom, it does not vanish at order in the IIB matrix model since S'efr[X] is 
shown to be nontrivial. Through the perturbative evaluation of S'efr[X] to all orders, we 
prove the finiteness of the X integration in eq. ( p^.30D for finite A^. 



2.3.1 One loop evaluation of 5'cfr[X] 

We substitute the one loop effective action SIq^°°^[X,C,] of eq. (|2.19|) into eq. ( p.30| ) and 
consider what kind of terms survive after ^ integrations: 

/«d{e--.V~'M = /«<^nil + ^ + (i(^)^ + ^)l(2-31) 

Here the products are taken over all possible different pairs of color indices (ij). When we 
expand the multi-products, we select one of the three different factors 1 or tr{S'^^j-^) /A or 
(tr{S^-j-^)/8 + (tr{S^-j'j))'^ /32) for each pair of (ij). Since the last two factors are functions 
of (x* — x^), they can be visualized by bonds that connect the "space-time points" and 
. More precisely, in order to remind us that the factors tr(S'^j^.^)/4 and (tr{S^-j^)/8 + 
(tr{S'll-j'j)y / 32) contain 8 and 16 spinor components of — we draw 8 or 16 bonds 
between and x^ depending on whether we take tr{S^-j-^)/4: or (tr{S^^j-^)/8 + {tr{S^-j^)y / 32). 
That is, each bond corresponds to each component of a spinor — We call these sets of 
8 and 16 bonds "8-fold bond" and "16-fold bond" respectively. In this way we can associate 
each term in the expansion of multi-products in eq. ( |2.31D with a graph connecting the space- 
time points by 8-fold bonds or 16-fold bonds. We don't assign any bond to the factor 1. 
Therefore the multi-products in eq. (|2.31| ) can be replaced by a summation over all possible 



graphs. 

Let us consider spinor component by component in order to discuss what kind of 
graphs survive after ^ integrations. Out of the 16 spinor components of we may focus on 
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a particular spinor component such as the first component ^J. We rewrite eq. ( ^.311 ) as 

„ .N-l N-1 

jdxj n ■ ■ ■ ^6 n c n(Co + c, ■ - iD). (2.32) 

i=l i=l i<j 

Here Co and Ci are functions of x^— and the other spinor components, ^2~^2; ' ' ' ; ^i6~^i6- 
We can then expand the multi-products regarding the factor Ci ■ — ^i) as a single bond. 
Then it is clear from the following considerations that the integrations of for all color 
indices i generate subgraphs associated with the first component of spinors, which connect 
all the N points without a loop. Such type of graphs are called maximal trees. 

i) Since has — 1 independent color components, the subgraphs having — 1 bonds 
remain after ^,1 integration. 

ii) If there is a loop in the subgraph, the contribution vanishes since a product of delta 
functions of grassmann variables on the loop vanishes: 

m - ema' - ^t) ■ ■ ■ - anm - er) = o. (2.33) 

We also note that all maximal trees contribute equally as we can see by performing 
integrations form the end points of the maximal trees. 

Therefore the whole integration of the fermionic degrees of freedom generates graphs 
which are superpositions of 16 maximal trees. However these 16 maximal trees are not 
independent. At each bond they are constrained to be bunched into the set of 8 or 16. From 
these considerations, we find a graphical representation of the one loop effective action: 

G:graph (ij):bond of G 

[(— ^) or {-{—f^f + —f^)]. (2.34) 



Here we sum over all possible graphs consisting of 8-fold and 16-fold bonds which can be 
expressed as superpositions of 16 maximal trees. For each bond {ij) of G we assign the first 
or the second factor depending whether it is 8-fold or 16-fold bond. 

Since each 16-fold bond contains 16 spinors, — S,^) and (x* — x^) must form Lorentz 
singlets by themselves, 

.l.M%))2 ^ ^(16)... _ ,n \ . (2.35) 

4^4^ 8 ^ ^ ^ (x^-x^)24- ^^-""^^ 

In the 8-fold bond, however, — ^•') and (x* — x^) couple as 

tr(5j^.)) = T^'^^Pix' - x^)^(x' - x^),(x' - x^)a(x' - x^)p/{x' - x^f\ (2.36) 
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Here 

[(e - ^')T'\{e - - ^^■)r'''.(e - e)] (2.37) 

is a totally symmetric traceless tensor as is shown in the appendix |^. 

If there were only 16-fold bonds, considerable simplifications take place since the graphs 
are reduced to maximal trees with 16-fold degeneracy and the interactions between the points 
only depend on their distances. We may cite the four dimensional model in eq. ( |2.21| ) as 
a simple model which shares such a property. The effective action of the four dimensional 
model can be represented as follows: 

dXd^ E n (2-38) 

Gigraph (jj):bond of G 

where we sum over all graphs whose bonds form maximal trees. Since ^ has four spinor 
components in four dimensions, 

(2.39) 

Therefore, in the four- dimensional model, the distribution of X becomes "branched polymer" 
type: 

Gimaximal tree (ij):bond of G ^ ' 

Note that all points are connected by the bonds, and each x*-' integration can be performed 
independently and converges for large x*-' on each bond where x*-' = x^ — xK Thus this system 
is infrared convergent. Although the integrations seem to be divergent for short distances, it 
is due to the failure of the one loop approximation since the theory itself is manifestly finite 
at short distances. We will consider the short distance behavior of the theory in subsection 
T^ . As we see in the appendix |C|, the dynamics of branched polymer is well known and its 



Hausdorff dimension is four. So this model constitutes an example of models for dynamical 
generation of space-time which predicts four dimensional space-time. As is shown in the 
subsequent discussions, the IIB matrix model is much more complicated. Nevertheless we 
expect that the structure of space-time is also determined dynamically in the IIB matrix 
model. 
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2.3.2 Infrared convergence 



In what follows we show that the X integral is convergent in the infrared region for finite 
to all orders of perturbation theory: 



J rfXe"^^«[^' < oo (infrared). 



(2.41) 



It shows that all points are gathered as a single bunch and hence space-time is inseparable. 
Indeed we can show that the multiple integral, eq. ( |2.41|) , is absolutely convergent. 

We use the following theorem to prove the above statement: An M -dimensional multiple 
integral is absolutely convergent when the superficial degrees of divergences of sub-integrals on 
any m- dimensional hyper-planes within R*^ are negative. Let us apply the identical linear 
transformations T to the variables X and ^ as 



(2.42) 



and consider the case where the values of ■ ■ ■ y"^ become large while y^+^ ■ ■ ■ y^ ^ are 



( 


] 










( 


] 




i 


\ 






= T 












= T 








) 






) 






) 






J 



fixed. We observe from the Feynman rules in subsection |2.1| that we can associate the factor 



(x' — x^ ] 



2 with every by absorbing half of the powers of the bosonic and the fermionic 



propagators which are connected to ^* — in eq. ( p.lOp . Although X dependencies of the 
diagrams are not entirely absorbed by this procedure at higher orders, the remaining factors 
make the infrared convergence properties even better. We explain this point by concrete 
examples at the end of this subsection. 

3 

We can reexpress (x* — x-')~2 and — in terms of y and rj through the identical linear 
transformations as {Y.kCky'')^^ aiid Y.kC^v'' respectively. The perturbative expansion of 
the effective action consists of many terms which contain fermionic variables through the 
combinations of ■ Let us consider the contribution from a term which contains 77^. If it is 
associated with the bond (ij), it means that Cl^ 7^ and we obtain the factor {Ci'y^ + ■ ■ ■)~2 
after the integration of r/J. We obtain the analogous factors from the integrations over all 
spinor components of fermionic variables up to 77™. Thus we can conclude that 

1 



^-Scff [?/] 



< 



rl6 pi 



(2.43) 



where are homogeneous functions of degree 3/2 of • ■ -y™. Hence the superficial degree 
of divergence is negative on this lOm-dimensional hyper-plane. Therefore the superficial 
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degrees of divergences are negative on any hyper-planes in R^*^*^^^^). This completes the 
proof that the theory is infrared finite for finite to all orders. Since the theory is manifestly 
finite at short distances, it also establishes that the IIB matrix model is a finite theory for 
finite N. 



Finally we give concrete examples which illustrate the fact that the factor (x* 



-3/2 



can 



be associated with every ^"^ — ^^ . At the one-loop level, the diagrams are made of the bosonic 
propagators ( p^.8| ), the fermionic propagators ( |2.9D and the ^ insertion vertex ( [^.lUI ). The 
Feynman diagrams contain them consecutively along the loop as depicted in the following 
figure: 




(2.44) 



Therefore we can diagrammatically understand that we can precisely assign the factor (x* — 



to {^'^ — in the one loop effective action. 



Since the expansion parameter is over the fourth power of the average distances be- 
tween the points, the infrared convergence is expected to become better in higher orders. In 
fact at higher loop level, there remain some vertices with nontrivial X dependencies which 
make the infrared convergence property even better after absorbing half of the powers of the 
bosonic and the fermionic propagators connected to every ^. For example. 




^2.45) 

is proportional to {C'')'^ /{^^^)^ ' {i'^^Y / {x'^^Y ■ (7^/((a;*-^)^(x*^)^), where the third factor is 
associated with the four gluon vertex at the center of this diagram and indeed makes the 
infrared convergence property better. Here we have introduced the notation = — 

anda^ = a-a- 



2.4 Short distances 

Until now we have considered the effective action for the diagonal elements which is valid 
when all of the diagonal elements are well separated from one another. It is analogous to 
the most generic moduli space in the super Yang-Mills theory although precisely speaking 
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there is no moduli in the IIB matrix model. It is natural to consider the second most generic 
moduli space next. Let us suppose that a pair of the bosonic coordinates are degenerate but 
the rest of the coordinates are well separated from one another and from the center of mass 
coordinates of the pair. 

We then decompose the N x N matrix valued variables into the 2x2 submatrices and 
the rest. The rest consists of the Hermitian (A^ — 2) x (A^ — 2) matrices and the complex 
2 X (A^ — 2) matrices. Let us further decompose the 2x2 submatrices into the center of 
mass part and the remaining SU (2) degrees of freedom. Since the off-diagonal elements of 
(A^ — 2) X (A^ — 2) matrices and the complex 2 x (A^ — 2) matrices are massive, we can 
integrate them out as before. We write the remaining degrees of freedom as 

\ 



X, 



(2.46) 



V 



e ) 



(2.47) 



,3 



N 



Ni 



(2.48) 



The resulting effective action can be written as 

S.,[4^^'\^P'^^'hx„t,xl. 

The first term is the SU{2) part of the original action. The second term is essentially 
identical to the SU{N — 1) case we have discussed in this section. The center of mass 
coordinates of the degenerate pair play the (A^ — l)-th coordinates. The only difference is 
that the strength of the potential between the center of mass coordinates and the others 
has doubled. The third term denotes the interaction between SU (2) part and the diagonal 
elements ofA^ — 2xA^ — 2 matrix. We neglect it since it is small compared to the first term 
in many cases as is discussed in the concluding section. 

We still need to consider the SU{2) part to determine the dynamics of the relative 
coordinates of the pair of the points. Here we can cite the exact solution for the SU (2) case. 
As we show in the appendix 0, the distribution for the relative coordinates is 
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We conclude that there is a pairwise repulsive potential of —8 Inr type when two coordinates 
are close to each other. It is clear that these considerations are valid for arbitrary numbers of 
degenerate pairs although the center of mass coordinates should be well separated. Although 
it is possible to repeat these considerations to the cases with higher degeneracy, the analysis 
becomes more complicated. Therefore we choose to adopt a phenomenological approach and 
assume the existence of the hardcore repulsive potential of the following form: 

S,,,,[X] = Y.g{x'-x^), (2.50) 

i<j 

where 

gix^ - x^) = I 7^^"^^^' - ""'^"'S) for (xj - x^ f « g 

\ for (x* - xJf > 5( ^ ^ 

In the following section, we investigate the structure of space-time by using the one loop 

effective action for the diagonal elements plus the phenomenological hardcore potential 5'core- 

Although our effective action is not exact at short distances, it presumably captures the 

essential feature of the IIB matrix model namely the incompressibility. If so our effective 

action may be in the same universality class of the full IIB matrix model. 

3 Possible scenarios of dynamical generation of space- 
time 

So far we have derived an effective action for the diagonal elements of A^^ and tp. We have 
identified the eigenvalues of with the space-time coordinates through the semiclassical 
correspondence as is explained in the introduction. Let us suppose that the eigenvalue 
distribution of A^ is extensive in d dimensions but supported only by finite intervals in the 
remaining 10 — d dimensions. Such a distribution may be interpreted that space-time is 
extensive in d dimensions but has shrunk in the remaining directions. The dimensionality 
of such space-time is obviously d. If such a distribution is derived from our effective action, 
we may conclude that the dimensionality of space-time is d. 

In the perturbative formulation of string theory, we need to assume a particular back- 
ground on which string propagates. Although certain consistency conditions are required 
for the backgrounds, we have no principle to pick a particular background from the others. 
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However we can in principle determine the possible structure of space-time from the IIB 
matrix model uniquely. Our effective action provides us an exciting possibility to realize it. 

In this section we study the effective action in order to see whether the IIB matrix model 
can explain the structure of space-time and its dimensionality in particular. As is explained 
shortly in more detail, our effective action is still rather complicated. One possible way to 
avoid it is to consider further simplified models which are supposed to contain the essential 
features of our effective action. Such models may contain several parameters which may not 
be easily determined. We therefore draw the phase diagrams of those models for all possible 
parameter regions and identify the structures of space-time in the respective phases. 

The aim here is to show the existence of models which can realize four dimensional 
space-time. We then need to show that the IIB matrix model indeed belongs to the same 
universality class. Although these investigations are still in progress, we find it very plausible 
that the IIB matrix model realizes four dimensional flat space-time. We explicitly propose 
some of possible mechanisms to realize realistic space-time. 

3.1 Models 

We consider the following ensemble in ten dimensions which is controlled by the 1-loop 
effective action ^] (|2.19|) and the core potential (|2.50|) : 



dXe-'^'-di J: n l(^) or (i(^y^ + ^)l. (3.1) 

Gigraph (ij):bond of G 

As is explained in section ^ we call the bond in association with tr{S(ij)Y ^^"^ ^^e bond 
with {tr{S(ij)YY /'i2 + tr{S{.ij)Y a 8-fold bond and a 16-fold bond respectively. The former 
gives 8 fermionic variables per bond while the latter gives a product of all 16 components 
and is proportional to {x'''^)~'^'^5^^^\^'^ ~ i^)- 

If only 16-fold bonds appear in a graph, ^ integration can be easily performed which gives 
a branched polymer as we saw in the previous section. The integrations can be performed 
bond by bond and the resultant interactions between a;*'s are functions of distances only, 
that is, (x* — for each bond. Thus we can estimate the partition function as 

dXe-^-rfee-"-"''[^'«] / dXe-'- ^ H 7^1^^' (3-2) 

G:maximal tree (ij):bond of G ^ ' 

where summation of G is over all maximal tree graphs. 
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When 8-fold bonds appear in the graph, the simple branched polymer picture breaks 
down. As is seen from eqs. ( p.36|) and (|2.37|) , tr(S'(jj))^ can be rewritten as 



tr(5(.,))^ cx C'^-Ap«....«s^^^ . . . s:iVi^^J{x^^Y\ (3.3) 

where ^"2 = Ca ~ = ~ (7M'^V"i-- a8 jg g^n invariant tensor. V^^^^p is a fourth 

rank symmetric traceless tensor constructed from x*-' : 

........ (x^^\^ .... .... 

T/'.? — rf^J rf^J rf.'^] rf^J _ ^ ' ( r^tJ ry-^.l _|_ -|- . . . 

^ {D + 2){D + 4) + + ---)^ (3-4) 

where D = 10. Therefore at each 8-fold bond there are many choices to select 8 fermionic 
components out of 16 and each choice gives different dependence which is specified by 
(jpuXpa^...o,>i_ Hence integration over ^ generally gives complicated interactions which involve 
many bonds and depend on relative directions of relevant bonds. An important point here 
is that the tensor V^^^y^p is traceless and vanishes if we naively take the orientation average 
for x^. Therefore it is natural to expect that the appearance of 8- fold bond is suppressed to 
some extent if the eigenvalues of are distributed uniformly in ten dimensions. 

As a simpler model which may capture the essential feature of our effective action, we 
may consider a model where we average the orientation dependence of the 8-fold bonds. 
Instead of considering the graphs consisting of 8-fold and 16-fold bonds, we introduce two 
independent maximal trees each of which maximally connects the N points with coordinate 
x^. When the two trees share a bond (ij), we regard it as a 16-fold bond and give the 
corresponding factor ((x*-')^ + g)"^"^- For a bond which is contained in one tree only, we 
regard it as a 8-fold bond and assign a factor ((x*-^)^ -|- g)~^ which represents the residual 
effect of the 8-fold bond after orientation average. Here we have introduced a short distance 
cut-off of order g in order to avoid the meaningless singularity at x*-' = 0. We also assume 
the hardcore repulsive potentials as given by 5'core at short distances: 

G;two maximal trees (jj):bond of G 

X [((x^^y + gr'e-^^' or ((x^^f + g)-'']. (3.5) 

Here the meaning of the factor e"'^/^ is as follows. As we have discussed in the previous 
paragraph, the orientation average suppresses the appearance of 8-fold bonds. Therefore we 
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can take the effect of angle-dependence at least partially by introducing a suppression factor 
e~^/^ for each 8-fold bond. We call this model a "double tree model" . 

If the IIB matrix model is really approximated by this model, A should be in principle 
fixed because the IIB matrix model has no free parameter. However to determine A precisely 
may be as difficult as solving the IIB matrix model itself. Therefore we treat A as a parameter, 
and examine the phase structure as A is varied form — oo to oo. We have employed numerical 
simulations. We show in the next subsection that two apparent phases exist in two 
opposite limiting cases A ^ oo and X ^ —oo . 

An important factor which is not taken into account in the double tree model is the fact 
that A is a sort of mean field and it may strongly depend on the state of the system, that is, 
the distribution of eigenvalues. Actually we will see in subsection |3.3| that A decreases when 
the eigenvalues distribute in lower dimensional space-time, and this effect favors to generate 
lower dimensional space-time. By introducing a rough approximation we can show that the 
free energy is minimized if the space-time is compactified to four dimensions. It appears 
that more detailed analysis of this effect is indispensable to arriving at our goal. 

3.2 Branched polymer phase and droplet phase 

In this subsection we analyze the double tree model defined by eg. ( [3.51) . It contains a 
parameter A and we consider two limiting situations. First we take A large enough. Then 
most of the bonds on one tree are bound to bonds on the other tree and the system behaves 
as an ordinary branched polymer made of 16-fold bonds only. We call this state a branched 
polymer (BP) state. As is reviewed in appendix 0, its partition function is given by 

A(0) = / d^x^j, (3.6) 

where we have introduced free energy per bond defined by Fbp = — ln(Q;c/i(0)). If we 
neglect the effect of intersections due to the core potential, we can set = e. We have also 
multiplied it by A^! since we distinguish A^ points x*. Xc is a short distance cut-off determined 
by the core potential. Since the Hausdorff dimension of BP is four, we cannot neglect the 
effect of intersections in dimensions less than eight and eq. ( |3.6| ) gives an overestimation of 
the partition function. In other words, Uc is reduced from e. Furthermore in dimensions 
less than four most of BP's cannot be packed and the number of possible graphs drastically 
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decreases. Without any other interactions besides 16-fold bond interactions, BP favorably 
expands in ten dimensions and we cannot get lower dimensional distribution of eigenvalues. 

On the other hand if we take A small, the number of 8-fold bonds increases and 16-fold 
bonds totally disappear in the A ^ — oo limit. Actually even in some positive values of 
A most of 16-fold bonds disappear, since the entropy increase caused by resolving 16-fold 
bonds into 8-fold bonds overcompensates the effect of positive values of A which favors 16- 
fold bonds. Our numerical simulation supports this fact as we will see later. Then the 
entropy is gained for such a configurations that as many as possible points gather around 
one another so that rearrangements of bonds occur easily. Indeed, without core potential all 
points are condensed into a finite area whose size does not depend on A^, as we show at the 
end of this subsection. If we take the core potential into account, they cannot condense into 
an infinitely dense state. Instead all of the cores are packed in an area whose size grows with 
A^. We call this state a droplet state. When a droplet state is realized in d dimensions, the 
size of the droplet is given by i? = 1{NY^^, where / is the core size. Since bonds flip locally 
in a droplet, we regard the droplet as a collection of clusters in each of which bonds can 
freely flip. The cluster is assumed to have n points inside. Then the volume of the cluster 
is f = I'^n, and the typical length of bonds in the cluster can be identified with the radius 
of the cluster r = In^^'^. Therefore each cluster is expected to contribute to the partition 
function by 

^(n) = (^ri!(^^)")%"e-""«/2, (3.7) 

where rig is the number of 8-fold bonds in the cluster. The first factor is the partition function 
of BP consisting of n points, and is squared since we sum over two independent maximal 
trees. We can then roughly estimate the partition function of the double tree model in the 
droplet phase by 

^droplet ~ {z{n)f/-^^ 



where Ag is the total number of the 8-fold bonds. Since d < 10, ^droplet is maximized by a 
small value of n which should be independent of A^. Therefore 2'dropiet can be written as 



^droplet 



Comparing the partition function of BP eq.( p.6| ) and that of droplet eq.( |3.9| ), we can 
tell that there is a phase transition or a cross-over between the two phases at some positive 
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Figure 1: Numerical simulation for N = 1000 . Horizontal axis is A. Stars indicate root 
mean squares of the eigenvalues and dots indicate the ratio of 16-fold bonds. For large values 
of A, the system is in a BP phase. For small values of A, it is in a droplet phase. 

value of A. For larger values of A, the system is in a BP phase and for smaller values of 
A, it is in a droplet phase. Our numerical simulation supports this transition (see fig. |l]). 
The nature of the transition is not clear at the present stage, although it has a similarity to 
the liquid-gas transition since we have no order parameter associated with a symmetry. The 
points distribute in full ten dimensions both in these two phases. Other possible interactions 
such as angle dependent interactions will modify these phases and we may obtain a phase 
with lower dimensionality. 

Finally in this subsection we see what happens if we do not have the core potential. 
In this case we can use mean field approximation to determine one-particle density in the 
following way. Let us denote p„ as the probability that a point is connected to n other 
points. Then the self consistent equation for one-particle density is given by 



where f{x) is the Boltzmann factor for each bond and given by (x^ + g)~^ in our case. Since 
N does not appear in these equations, the one-particle density does not depend on A^. Hence 
all points are condensed into a finite area V which does not depend on A^. We call this state 
a mean-field state. We can find a lower bound of the partition function by replacing the 
length of each bond with the size of the whole system: 




(3.10) 
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Since eq. ( ^.111 ) contains extra A^! compared with Zbp in eq. ( |3.(j| ), we have to increase A as 



log — i> oo in order to keep the state in BP phase. Thus we find that any state falls into a 
mean-field phase if we do not have the core potential. 

3.3 A mechanism favoring lower dimensions 

Although the double tree model has taught us a lot about the IIB matrix model, we may 
need to extend it in order to be more faithful to our effective action. As we saw in subsection 



3TT| , the interaction of 8-fold bonds generally depends on the relative angles of the vectors . 
If eigenvalues extend in ten dimensions, averaging over directions gives a suppression factor 
to each 8-fold bond, because the tensor V^i^xp is traceless. This suppression becomes weak 
if eigenvalues collapse into lower dimensions. Therefore if there are considerable number of 
8-fold bonds, the partition function is suppressed in higher dimensions. However the entropy 
which is associated with possible orientations of bonds favors higher dimensions. Therefore 
there is a competition between the two and it is conceivable that four dimensional space- 
time is realized as a compromise. In what follows we roughly estimate the partition function 
assuming that the system is in a droplet phase and show the possibility of this scenario. 
Assuming that lie isotropically in a d dimensional subspace, we can estimate the 

ij\4: 



orientation average of V^l^ as 



-{5p.u5\p + S^xSup + ■ ■ ■) 



d{D + 4) 

+ p^:^yp^:^('^Mi''^Ap + 5MA<^^p + ---)- (3-12) 

Here S stands for the projection operator to the d dimensional space-time while S is the 
Kronecker's delta in the original D = 10 dimensions. The interactions of 8-fold bonds in 
general involve many V^j^Ap's whose Lorentz indices are contracted in various ways. It is not 
easy to estimate the combinatorics which come out from such contractions exactly even after 
such simplification as eq. (|3.12|) . It may be qualitatively valid to replace 6 by 6. Then the 
relevant combinatorics become independent of d. Under such an approximation, each 8-fold 
bond becomes to have a wight effectively which is proportional to 

e-^/^ - f^M = (p + 2)P + 4) - id + 2)id + A) ^lidh- dTi)) ' ^^''^^ 
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Figure 2: Free energy as a function of dimensions. We take 7 = 17/6. Free energy is 
minimized at four dimensions. 

The weight fes{d) is a decreasing function and favors lower dimensions. Competition 
between this weight and entropy which favors higher dimensions determines space-time di- 
mensions in this scenario. If we assume that all bonds are 8-fold and the entropy per bond 
increases as d'^, the partition function for a droplet can be expressed as 

^droplet = iV!e-2^^('^) ~ iV!(/eff(rf))2^d2^^. (3.14) 

The free energy F has a minimum at ci = 4 when 7 is between 2.6 and 3.1. Fig. ^ shows F 
as a function of d. If this is the mechanism to generate our space-time, we need to determine 
detailed numerical factors to prove it. We stress here that the model is already fixed by 
nature and we don't have freedom to tune parameters. 

3.4 Other mechanisms for understanding space-time structure 

There are two other possible mechanisms to generate lower dimensional structure of the 
space-time. 

One way is to approach from a BP phase in which we regard 8-fold bonds in eq. (|3.1| ) 
as perturbations to BP's consisting of 16-fold bonds only. Since Hausdorff dimension of 
BP is four, there is a possibility that small perturbations might compress the system into 
four dimensions. If we neglect the effect of self-intersections and naively confine BP into d 
dimensional space, mean density becomes zero for d > 4 and diverges for c? < 4 in large 
limit. Then for d < 4, the core potential extends the bond length and energy drastically 
increases. In other words most of the topologies of BP are prohibited to avoid the increase 
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of energy and the entropy is decreased very much. For d > 4, large distances between the 
points prevent 16- fold bonds from resolving into pairs of 8-fold bonds, and we cannot gain 
the entropy of having 8-fold bonds at various places on the polymer. On the other hand at 
d = 4, 8-fold bonds can move around on the polymer and the system behaves as a gas of 
8-fold bonds. Therefore only at (i = 4 we can gain the entropy of 8-fold bond gas without too 
much excess of energy due to the core potential. The mechanism discussed here is based on 
the picture of BP and seems totally different from the one discussed in subsection ^.3| which 
is based on the picture of droplet. In four dimensions, however, these two pictures can be 
complementary, since the size of space-time predicted in these pictures are both N^^^l and 
agree with each other. 

Finally we point out a possibility that four- dimensional space-time is realized in the in- 
termediate region of fig. ^ even for the naive double tree model where dimension dependence 
of A is not taken into account. Here we consider dimension dependence of the ratio of 8-fold 
bonds, rg. In subsection we saw the transition between two limiting phases, namely a BP 
phase for large A and a droplet phase for small A. For intermediate values of A, the partition 
function can be expressed as 

Z(A) ~ AT! e-^(^(^)+^'-«), (3.15) 

where F{\ = oo) = Fbp and F(A = — oo) = -fdropiet- Obviously the free energy F{X) defined 
above is always smaller in higher dimensions. Therefore the system always expands in full 
dimensions, if rg is independent of dimensions. Our numerical simulations, however, showed 
that rg is smaller in lower dimensions for fixed values of A. It is consistent with the intuition 
that each point has a smaller number of neighbors in lower dimensions and the probability 
that two 8-fold bonds coincide to form 16-fold bond is larger. If the dimension dependence of 
Arg exceeds that of F{X), there is a possibility that lower dimensional phase becomes stabler 
than the ten dimensional phase for some value of A. 

4 Conclusions and Discussions 

In this concluding section, we first estimate the magnitude of the higher loop contributions 
to the effective action. Since our expansion parameter is g"^ over the fourth power of the 
average distances between the points, long distance behavior is well described by the one-loop 
effective action 5*^^ '""''[X]. However there is a possibility that small contributions may pile 
up at each loop since color indices circulate around the loop. Our strategy in this section 



24 



is to assume that the one loop level effective action determines the distribution p{x) of the 
eigenvalues. We then use p{x) to estimate the magnitude of the higher loop effects. This 
procedure is self-consistent if we find that the higher loop effects are small. Of course, our 
perturbative expansion breaks down if the expansion parameter grows as a positive power 
of N. What we show here is that the higher loop corrections are indeed finite in the large 
limit for finite coupling constant, and our approximation based on the loop expansion can 
be fully justified. 

Although we have argued that higher loop effects improve the infrared convergence prop- 
erties in general, there are diagrams which possess the identical infrared convergence prop- 
erties with the one loop diagrams in our estimation procedures just explained. For example 
we may consider the one loop self energy corrections to the propagators of the off-diagonal 
components. We then consider the insertions of them into the propagators of the one loop 
diagrams which contribute to the effective action. It is a part of the two loop level contri- 
butions. 

Recall that the matrix model (|1.1| ) has been constructed by the large N reduction pro- 
cedure from the SU{N) gauge theory. Thus, in the leading order in expansion, the 
perturbative expansions in g"^ in the matrix model are identical to those in the gauge theory 
were it not for the fermionic zeromodes. They are precisely the part of the quantum correc- 
tions we quoted above. Then the ratio of (/ + l)-loop contributions to /-loop contributions 
may be estimated by dimensional analysis as 

k 



^ Sd^^xp{x) ' ^ ' 

where p{x) = Y,k ^{^ ~ ^'') is a single-particle distribution of the eigenvalues. 

We first assume that x is distributed uniformly in a d dimensional manifold with size R. 
From the correspondence with gauge theories, we can say that there is no divergences from 
small x regions in ci > 4, and from large x regions in c? < 4. Thus, 

"^^d -4_ / R'''^ d>4: 
a X X — < ^^_4 ^ ^ ^ 5 l^"^J 
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where / is a typical size of short distance cutoff. Therefore, 

where we used the fact that the core size is I ~ g^^"^- Since ~ {R/lY, we may conclude 
that the quantum corrections are divergent ioi d > A but finite ioi d < 4 in the large 
limit. These arguments seems to imply the logarithmic divergence for d = A. However here 
we need to take into account of supersymmetry. It is well-known that dimensional reduction 
from ten dimensional super Yang-Mills theory down to four implies J\f = 4 supersymmetry 
in four dimensions. Furthermore four dimensional super Yang-Mills theory with A/" = 4 
supersymmetry is also known to be finite. So we expect that the quantum corrections for 
d = 4 case is finite. 

We have found a divergence in the above argument for d > A. However the situation 
is totally different if we assume that space-time is branched polymer like for d > 4. We 
recall that R ~ N^/^l for such space-time since its Hausdorff dimension is four. If so, the 



estimation of eq. ( [4.1[ ) is essentially the same as the four dimensional case, and we may 
conclude that quantum corrections are finite also for d > 4. Thus we can show that the 
corrections coming from planner diagrams are finite in the large N limit if the coupling 
is kept fixed. 

Next let us consider higher order effects of expansion. It is clear that the leading 
l/N"^ corrections contain two less summations over color indices compared with the planar 
contributions. However they do contribute since they are at least finite. In fact they can be 
of the same orders of magnitude with the planar contributions since they are more singular at 
small X regions. For example let us consider a three loop correction to the gluon self-energy. 
Such a nonplanar contribution should go like 

q^y i i ^ . (4.4) 

(x* - x^Y {x^ - x^Y {x^ - x^Y ^ ' 

Here the important point is that the superficial degree of divergence for the integration vari- 
able x^ is reduced compared with that of the corresponding planar diagram. For dimensions 
up to four, the above expression may be evaluated as follows: 

/ d^^xp{x) 

= 9'NjJj^„ (4.5) 
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where we set (x* — x^) = I ~ g^^'^ to estimate the magnitude of the correction at short 
distances. Since ~ {R/lY, the contribution is again nonvanishing and finite in the large 
N hmit. 

Therefore infinite numbers of diagrams contribute in dimensions up to four. This means 
that all orders of expansion contribute equally. On the other hand the nonplanar 
contributions diverge in general for d > A although they are much smaller compared with 
the planar contributions unless space-time is branched polymer like. If the space-time is 
branched polymer like, the situation is essentially identical to the four dimensional case 
since its fractal dimension is also four. So we may conclude that all orders of expansion 
contribute in all dimensions. 

Our findings may be interpreted that g^t ~ 1 since we have argued that string coupling 
constant is inversely proportional to . Apparently the double scaling limit is naturally 
taken in the IIB matrix model. However our investigations here also imply that the string 
scale {a'Y ~ 5'^- It is in disagreement with our previous estimate {a'Y ~ q'^^ based on 
the loop equations. We believe that reexaminations of the loop equations should be able to 
reconcile this discrepancy. 

We also would like to comment on the universality of IIB matrix models. It is clear from 
the investigations in this section that the renormalizability and finiteness of four dimensional 
super Yang-Mills theory is crucial to control the quantum corrections in the IIB matrix 
model. If we were to add higher dimensional terms to the IIB matrix model, we may need 
to renormalize the theory nonperturbatively. We expect that such a renormalized theory 
belongs to the same universality class with ours. Therefore we may argue that the IIB matrix 
model is universal and our Lagrangian may correspond to the fixed point Lagrangian. 

In this paper, we have derived an effective action for the IIB matrix model which is 
valid when the eigenvalues of are widely separated. It is the effective action in terms of 
the super coordinates of space-time. It has clearly shown that the theory has no infrared 
divergences and the universe never disintegrates. The exciting possibility is that we can 
determine the dimensionality of space-time by studying such an effective action. 

There are order A^ pairwise attractive potentials between the space-time coordinates. 
Such interactions may be classified into the 16-fold bond and the 8-fold bond types. If there 
were only the 16-fold bond interactions, the space-time coordinates form branched polymer 
type configurations with the fractal dimension four which stretch in ten dimensional space- 
time. Such a distribution of the eigenvalues cannot be considered as continuous space-time. 
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We have found that the closely packed distributions of the eigenvalues can be realized when 
we take the 8-fold bond interactions into account. Such distributions can be interpreted as 
continuous space-time. Here the existence of the hard core repulsive potentials between the 
eigenvalues is also very important. Although more detailed investigations are required to 
show that four dimensional space-time is realized in this model, we have proposed possible 
mechanisms to realize the realistic dimensionality. 

We believe that what we have already learned from the IIB matrix model is very sig- 
nificant. It is not defined on the fixed background metric but space-time is dynamically 
determined as the vacuum of this model. The extent and the dimensionality of space-time 
can be obtained from the eigenvalue distributions of A^. In this model, the vacuum is not 
empty and may be classically represented by the diagonal matrices. The matter (for ex- 
ample D-objects) arc the local fluctuations which float on the vacuum. In this way the 
space-time and the matter are inseparable and they determine each other. Note that such a 
picture never emerges as long as we deal with the effective Lagrangians for a finite number 
of D-objects. 

If four dimensional space-time is generated dynamically, it will be plausible that the 
space-time has four dimensional Lorentz invariance. Also it will be intrinsically flat. Here the 
core potential plays an important role to protect the space-time from shrinking. On the other 
hand long-distance dynamics protects it from expanding infinitely. This means vanishing of 
the cosmological constant. In our scenario, stability of generated space-time guarantees 
absence of cosmological constant in an effective theory of gravity, that is, Einstein gravity. 
If vacuum energy is generated, eigenvalues are rearranged so as to restore stability of space- 
time. It is dynamically tuned to be zero. The existence of the hardcore repulsive potential 
at short distances and the infrared finiteness are responsible to achieve this remarkable 
feat. It does not explicitly depend on the supersymmetry and hence should be independent 
of whether supersymmetry is spontaneously broken or not. We believe that this fact alone 
constitutes a major achievement and hence underscores the validity of the IIB matrix model. 

We recall the M — 2 SUSY transformations for the zeromodes: 

= 0, 

= ieir^e, (4.6) 

and 
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= 0. (4.7) 

The linear combinations of the above translations 6^ = 6^ ± 6"^ form the Af = 2 supersym- 
metry algebra. The vacuum expectation value of the supersymmetry transformation of ^ 
must vanish if supersymmetry is not broken by the vacuum. It is clear from eg. ( |4.7|) that 
the M = 2 supersymmetry of 5^ are broken completely. It is in contrast to the D-string 
case where the half of the supersymmetry can be preserved. If so, this model realizes the 
vanishing cosmological constant without supersymmetry as we have argued. 

It is also possible to write a scenario to obtain realistic gauge groups in this model. 
Although we have considered the most generic case where all eigenvalues of are different 
from each other, we may assume that some of the eigenvalues remain degenerate. Then the 
vacuum configurations are represented by the block diagonal matrices where each block is 
m dimensional. Then the gauge group SU{m) is realized (note that U{1) part is used to 
make the space-time coordinates). It is also possible to realize the standard model gauge 
group in an analogous way. We may again integrate out the off-diagonal elements and 
obtain the effective action for the U{m) submatrices. The resultant effective action must 
closely resemble the gauge theory with the Planck scale cutoff since we have the local gauge 
invariance as the manifest symmetry of the effective action. 

Since we have argued that it is possible to obtain the realistic gauge groups in four 
dimensions in this model, it is natural to ask to what kind of string compactification it 
corresponds. It presumably corresponds to the compactification into the interior of that 
is B^. It may be possible to play the analogous games with Calabi-Yau compactifications of 
heterotic string. It is possible to obtain chiral spinors in the fundamental representations of 
the gauge group in four dimensions from the chiral spinors in the adjoint representations in 
ten dimensions. For this purpose we need to assume that there is nontrivial gauge field con- 
figuration in which produces a nontrivial index for the Dirac operator. We immediately 
think of a two dimensional example such as magnetic vortices. Then the tensor products of 
three of them are the possible candidates for such nontrivial gauge configurations. 

We also note the similarity between our effective Lagrangian and the dynamical triangu- 
lation approach for quantum gravity. Due to the existence of the hardcore potential of the 
Planck scale, the unit cell of the effective action is also the Planck scale. The pairing interac- 
tion between the eigenvalues may be identified with the bonds in the dynamical triangulation 
approach. The integration over the fermionic zeromodes ensures that the model sums the 
contributions with all possible connectivities of the bonds just like the dynamical triangula- 
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tion approach. This analogy may be useful to reconstruct the metric out of the IIB matrix 
model. We also conjecture that the general coordinate invariance is present in this model 
due to the following reasoning. There is the permutation symmetry Sn which permutes the 
color indices. It is a subgroup of the full gauge group SU {N) and an exact symmetry of our 
effective action. Since it does not change the density of the eigenvalues, it should be part of 
the volume preserving diffeomorphism group in the continuum limit. Similar argument has 
been made in the dynamical triangulation approach. However it suffers from the infamous 
conformal mode instability while the IIB matrix model is free from such a disease which is 
again a remarkable merit of this model. So the gauge and the diffeomorphism invariances 
may be unified into the SU (N) symmetry of the IIB matrix model. 

Although the investigation of this model is still in the initial stage, this model turns out to 
be a finite theory which is free from both short distance and long distance divergences. The 
existence of such a theory itself is very impressive. It is a manifestly covariant formulation of 
superstring which enables us to determine the structure of the vacuum, namely space-time. 
We have argued that it can solve the problems such as the cosmological constant problem 
which appear to be insurmountable before the advent of it. Therefore we believe that it will 
reveal further truths concerning the structure of space-time and the matter. 
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A Fierz transformation 



In this appendix, we investigate the totally symmetric tensors which can be constructed out 
of a ten dimensional Majorana-Weyl spinor C,- The effective action contains such a term as 
trS^ where 5*^*^ = ^r'^"'^^ x^jx^ and n is an even integer up to 8. Therefore we are interested 
in totally symmetric tensors. 

We first point out that the only nonvanishing tensors which are quadratic in ^ are 
Thus the Fierz transformation is performed quite easily, and we obtain the following identi- 
ties: 

= 0, (A.l) 
From eq. ( [A.l|) and eq. ([A.2|) , we can prove the following identity: 

(er%0(er-^^0(er^V) = 0. (A.3) 

proof: Since the first factor of the left-hand side is antisymmetric in a and f3, the other 
factors can be antisymmetrized as the right-hand side of eq. ( |A.2D . Therefore using eq. (|A.2| ) 
we have 

L.H.s. = -i/2(eT'^"^0(er^\0(erV0, (a.4) 

and it vanishes after the contractions of a and f3 due to eq. (|A.1| ). (q.e.d.) 

Eq. ( |A.1| ) shows that any totally symmetric tensor made of four ^ vanishes. Since sixteen 
C, form Lorentz singlet, we can see from the "interchange of particles and holes" that any 
totally symmetric tensor made of twelve ^ must vanish. From eq. ( |A.3| ), we can also conclude 
the absence of the totally symmetric tensors made of six and ten spinors in an analogous 
way. 

Thus we have shown that the totally symmetric tensors made of ^ are exhausted by the 
following hst: {1, T^^^p = {^T^^'f^Oi^^'^'^'yOi^^^^sOi^^''^ aO^ ^^^}- It is easily checked that 
j^fj-i^Xp jg totally symmetric traceless tensor in the following way. Similar use of eqs. (|A.2| ) 
and ( |A.3|) to the above gives the antisymmetric part of T^^^p as 

= 0, (A.5) 

and from eq. ( |A.1| ) we immediately see that T^^\ = 0. Therefore T^'^^^ is totally symmetric 
and traceless. 
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B SU(2) Matrix Model 

In this appendix, we solve supersymmetric SU(2) matrix models in various dimensions, 
D = 3,4,6,10: Q 

where Af^{fi = 0,1, D — 1) and ip are 2x2 Hermite matrices, and ^p is a. spinor in D 
dimensional super Yang Mills theory. Each spinor consists of real components, where po 
is 2,4,8, 16 for D=3,4,6,10, respectively. We denote D-dimensional alpha matrices by a^'s 
{a'^ = r^r'^) and the first three of them can be represented as follows: 



= ® 

r I2 for i=0 

p'=la' for 1=1 . (B.2) 
I for i=2 

Performing the integral over fermionic variables gives the pfaffian 

Pf(a^A^) = Pfia^'T^A^^), (B.3) 

where T°- is the generator of SU{2) in the adjoint representation, namely {T°-)hc = ei,ac- Now 
we "rotate" by Lorentz transformation so that only Aq, Aiand A2 are nonvanishing. Eq. 
( p.3|) then reduces to the three-dimensional calculation, 

{Pf(p*TM»)f-/^ ~ {e^.keabcAtA'^Al}"^/', (B.4) 

which corresponds to (detA)^-^/^ when we regard A^ as a 3x3 matrix. 

Next step is the integration over three ten-dimensional vectors ^4°, which we reduce to 
the integration over three-dimensional vectors. The Jacobian for this reduction is the volume 
of the parallelepiped spanned by the three vectors, A^, A? and A? , to the [D — 3)-th, which 
is nothing but | det 

To estimate the behavior of the integral, we take the following parametrization for three 



vectors A\, 

t\ ( a\ ( h 



OJ'\yJ'\z 



(B.5) 



^Some results presented here has been aheady reported in Rcf. |g2[. Similar calculations have also been 
done in a different context in Ref. [E3|. 



^Although the integration does not lead to a pfaffian in D = 6, the result (B.4) holds in this case as well. 
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where y and z are two dimensional vectors while is two dimensional zero vector. After 
these considerations, we obtain 



dAdipe ^ 

= J drr^dadbd^yd\{r{y x z)y^^\r{y x z)\^-^ 

X exp{-^[r^{y^ + z^) + {az - hyf + {y x zf]). (B.6) 

Integrating over a and b yields 

I dadhe-^'^^-'y^" = {z^y^ - {z ■ yfy'^' = \y x z\-\ (B.7) 

Next, let us integrate over y and z. For D=3 case the integral becomes identically zero, 
since {r{y x z)} takes the both signs. For the remaining dimensions, we can estimate the 
integration when r is large and small respectively, which suffice the present purpose. When 
r is large, the ffist term of the argument of the exponential in eq. ( [B.6|) becomes dominant. 
We can perform the integrations over y and z by rescaling them as y/r and z/r. In this 
way eq. (p.6|) is estimated to be 



J r^dr r ^ for D=4 

Jr^drr-^^ for D=6 (r^ > (B.8) 

/ r^dr r~^^ for D=10 

When r is small, the ffist term of the argument of the exponential becomes negligible. We 
can perform the integrations over y and z independently from r. Eq. ( p.6|) is estimated in 
this case to be 



/ 



J r^dr for D=4 

Jr^drr^ for D=6 (r^ < (B.9) 

/ r^dr for D=10 

Finally we confirm that the long distance behavior of r integral (|B.8| ) agrees with the 
result of the perturbative calculation in the text. Recall that the one- loop effective action 
for r and is written as 

Sl^''"'y,^]=trHl + Sn, (B.IO) 

where 
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Thus the integration over ^ 

J rf^rrf^'^^e-^^ffl'^'^l (B.12) 

leads to eq. (^). 

C Branched polymer 

In this appendix we give a brief introduction to thermodynamics of branched polymers (BP) 
without self-avoiding effect in D dimensions, (see for example 0]). The partition function 
of BP with points and (A^ — 1) bonds is given by 

iV-l oo 

zn=t. Ud''^' n fi^'-^niitt, (C.l) 

BP i=l (ij):bond n=l 

where the summation is taken over all possible topologies of branched polymers. t„ denotes 
a positive weight assigned for the points to which n bonds are connected and denotes 
the number of such points in a given configuration. The weight function for bonds f{x) can 
be Fourier transformed as 

fip) = J d^xf{x)e-'^P- = /(0)(1 - C2P^ + C4/ + ■■■), (C.2) 

with a positive coefficient C2. Since we are interested in the thermodynamic limit where N 
goes to infinity, we consider a grand canonical partition function: 



Z=Y: NZ^ki^-'l (C.3) 

N=l 



In order to demonstrate the scaling behavior, we consider a two point correlation function. 
If we pick up a pair of points in a particular configuration, they are uniquely connected by n 
bonds. If we assume that they are separated by a fixed distance, n varies from a configuration 
to another. Therefore we can express a two point function in the momentum space as follows: 

G{p,ko) = |d^(x-i/)G(x-y,Me'^^""') 
00 J. 
oc Y.kr'iti + t2b + ^^b' + ■■■)'{ f{p)rB'^-\ (C.4) 

n=l ^• 

A graphical illustration of this equation is given in fig.^. The factor (ti + + ■ ■ ■) is the 
contribution from each end point in fig. ^ B represents each blob in fig. ^ except those 
at the end points, b is graphically represented in figfl. It satisfies the following Schwinger 
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Figure 3: Graphical representation of a two point function. A pair of points in a given 
configuration can be connected by bonds in a unique way. In this figure they are connected 
by five bonds. Such bonds form a random walk type object. The remaining points in a 
branched polymer are connected to the points in this object. They are represented by the 
blobs in this figure. 



b = 



Figure 4: Graphical representation of the Schwinger Dyson equation for b. White circles is 
not associated with any wight. 



Dyson equation 



b=kofmti + t2b+^^b^ + 



kh{b), 



where k = kof{0). The blob B can be related to h{b) as 



B = h + hb + ^b' + --- = h'{b). 
From these considerations, we can see that G{p, ko) can be written as 

Gip, ko) oc — - — . 

l-kofip)h'ib) 

We also note that Z'{b) = h{b), since the partition function ( |C3|) is given by 

Z{b) = 1 + ti6 + t26V2! + ta&VS! + ■ • •■ 

Since the averaged number of the points in the branched polymer is given by 

k dZ kh{b) 



N 



Zdk z{b)(f. 



(C.5) 



(C.6) 



(C.7) 



(C. 



(C.9) 



we must tune the fugacity k^ so that dk/db approaches zero from a positive value. Eq. 
solves as a function of 6: k = b/h(b). A typical solution is illustrated in fig.^. Generally 
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Figure 5: Typical relation between k and b is drawn. A black dot in the figure indicates the 
critical point. The large N limit is taken by approaching this point from the left. 



ti ^ and some of t„'s with n larger than 2 are non-zero. We observe that k{b) vanishes 
at 6 = and b = oo. Since it is positive definite, there is a critical fugacity kc where an 
averaged number becomes infinite. 

Near the critical point, k{b) can be approximated by 

k^k,-^ib-b,f, (C.IO) 

where c is a positive constant determined by the set of t„'s. The partition function is given 
by the integral of h{b) over b and it behaves near the critical point as 



Z = const. (C.ll) 

kc\ c 

The universal part (the second term) determines the power —3/2 in the following expansion: 



Z=J2 N-^'^ ( f ) ~ E iV-3/2e-^^. (C.12) 

This is because the universal part of the infinite sum of the r.h.s. in eq. ( |C.12| ) can be 
estimated by an integral / dNN~^^'^exp{—{kc — k)N/kc). Therefore partition function for a 
fixed is = N~^^'^(f(0)/kc)^. The leading non-universal behavior is given by 

^iv ~ (/(0)«c)^ (C.13) 

where = ^/kc takes a value of order one. In a special case where tn = 1 for all n, we find 
that k = be~^ and the critical values are kc = 1/e and be = I. This is the case relevant to 
our analysis. 

Inserting kh'{b) = 1 — h{b)k'{b) (which is derived from eq. ( |C5D ) into eq. ( |(J.7] ), we 
obtain a scaling behavior of a two point function near the critical point. For small p, it 
behaves as 

kLh'{bc)f{p) 



G{p, ko) oc 



(/(P)//(0))(1 - ^2cikc-k)hibc)) 
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where = \/2ch{hc) / C2- In order to take a scaling limit, we have introduced the physical 
momentum pphy = p{kc — k)^^^^. In this way, we have obtained the scaling relation: 

(fc^ _ k)-'/'G{x{k, - k)-^l\ ko) '-^^ g{x) (C.15) 

where g{x) is a certain function. Since the dominant contribution comes from ~ kc/{kc — 
k), we can observe from eq. ( |C.15|) that the Hausdorff dimension of a branched polymer is 
four. 
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